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ABSTRACT 

The generation and evolution of small amplitude long wavelength travelling distur- 
bances in rotating disk flow is the subject of this paper. The steady rotational speed of the 
disk is perturbed so as to introduce high frequency oscillations in the flow field. Secondly, 
we introduce surface imperfections on the disk such as roughness elements. The interac- 
tion of these two disturbances will generate the instability waves whose evolution is 
governed by parabolic partial differential equations which are solved numerically. It is 
found that, for the class of disturbances considered here (wavelength on the order of 
Reynolds number), eigensolutions exist which decay or grow algebraically in the radial 
direction. However, these solutions grow only for frequencies larger than 4.58 times the 
steady rotational speed of the disk. The computed receptivity coefficient shows that there is 
an optimum size of roughness for which these modes are excited the most. The width of 
these roughness elements in the radial direction is about .lifc where r 0 * is the radial 
location of the roughness. It is also found that the receptivity coefficient is larger for a 
negative spanwise wavenumber than for a positive one. Typical wave angles found for 
these disturbances are about -26°. 
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1. INTBQDUCTIQN 

In this paper, we investigate the receptivity and non-parallel stability problem for 
travelling disturbances in a rotating disk flow, A significant feature of the boundary- 
layer flow on the rotating disk is that its velocity component in the radial direction, the 
crossflow velocity, has an inflectional point. This causes an instability, often called 
crossflow instability, first noticed by Smith (1946) and investigated theoretically and 
experimentally by Gregory et al. (1955). When the crossflow component is combined in a 
particular direction with the velocity component in the azimuthal direction, they form a 
mean velocity profile which has an inflection point at which the velocity is zero. This 
permits an inviscid neutral disturbance with zero frequency which connects to unstable 
modes at other wavenumbers (see Stuart in Gregory, Stuart & Walker (1955)). The 
stationary disturbances appear as vortices, called crossflow vortices, which rotate in the 
same direction. This phenomena is also observed in other flow geometries such as 
rotating cones and swept wings, and it is this latter practical application which generates 
the interest in the problem since it is found that the crossflow instability dominates the 
boundary-layer transition process near the leading edge of a swept wing. 

In the rotating disk flow, the crossflow vortices spiral outward at an angle whose mag- 
nitude is approximately 10°-14°, with respect to the azimuthal direction (Fig. 1), and these 
vortices first start to appear at a Reynolds number of about 286. Wilkinson and Malik 
(1983), using hot-wire techniques, mapped out the complete wave pattern on the disk and 
found that the stationary disturbances originate from isolated roughness sites on the disk. 
Mack (1985) computed the wave pattern observed in the Wilkinson-Malik experiment 
using the linearized stability equations and assuming a white spectrum at the source of the 
wave pattern. Malik (1986) calculated the neutral curve for stationary disturbances, and 
found numerically another stationary viscous mode which corresponded to zero wall- 
shear stress of the mean velocity profiles. At high Reynolds numbers on the lower branch 
of the neutral curve, the wavenumber a behaves like a = const/B 1/2 , and at the upper 
inviscid branch they behave like a ~ const. Here a is the nondimensional wavenumber in 
the radial direction, which is nondimensionalized by the length scale (v/ f2 0 ) 1/2 , and R is 
the Reynolds number, defined by R = r*(£2 0 / v) 1/2 , where v is the kinematic viscosity, D 0 is 
the angular velocity of the disk, and r* is the radius where the local stability analysis was 
done. 

Hall (1986) developed an asymptotic analysis to investigate the two branches of the 
neutral curve. He showed that the inviscid mode is characterized by a two-layer structure. 
The upper deck is governed by the Rayleigh equation, and a thin layer near the wall is 
needed to satisfy the no-slip boundary condition at the wall. The lower branch neutral 
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curve is characterized by a triple-layer structure and has a long wave length. This mode 
has its structure fixed by a balance between viscous and Coriolis forces. At high Reynolds 
numbers, the lower branch neutral waves are inclined at about 39° to the radial direction. 
It is also observed that the formal asymptotic theory results, and the linear stability 
calculations based on parallel-flow assumption, agree for a large range of finite Reynolds 
number. 

In a recent paper by Balakumar & Malik (1990) (referred to as BM below), linear 
stability characteristics of stationary and travelling disturbances in a rotating disk flow 
were presented. Figure (2) shows the neutral stability curves in the radial wavenumber a 
vs Reynolds number R plane for different nondimensional frequencies. Here the nondi- 
mensional frequency o) is defined as the ratio of the dimensional frequency of the distur- 
bance (o * to the rotational speed of the disk D 0 . For a detailed discussion of this, the reader 
is referred to BM. In this paper, we are concerned with the lower branch neutral curve at 
high frequencies. It is seen from Fig. (2) that at high frequencies, (0 = 6.0 and 7.9, the 
wavenumber a varies like R 1 . The waves for these modes are inclined at negative 
angles, see Fig. (1). For a> = 7.9, this angle is about -35.32°. 

The stability results in BM have been obtained from linear stability theory with paral- 
lel flow assumptions being used. In the rotating disk case, this assumption amounts to 
locally replacing the variable radius r* by the constant radius r*. This approach is justi- 
fied if the wavelength is much smaller than the radius r*. This requirement is violated 
when the wavenumber is of the order of J? 1 , and the validity of the Orr-Sommerfeld type 
approach becomes questionable. This is similar to the situation which exists in the Gortler 
problem, Hall (1982,1983). It was shown by Hall that, at 0(1) wavenumbers, the Gortler 
instability is governed by a parabolic type partial differential system. In this paper we 
investigate this R 1 BM mode of instability at high frequencies and allow for nonparallel 
effects in a self-consistent manner. 

We formulate the problem as a receptivity problem similar to that analyzed for the 
Gortler case by Denier, Hall, & Seddougui (1990). We perturb the steady rotational speed of 
the disk so as to introduce high-frequency oscillations in the flow field. Secondly, we 
introduce imperfections on the surface as roughness elements. The interaction of these two 
flow effects, the unsteady one due to the oscillations of the disk and that introduced by the 
roughness, will generate the instability waves whose evolution is governed by parabolic 
partial differential equations. These partial differential equations are solved numeri- 
cally. In Section 2 we describe the formulation of the receptivity problem for the rotating 
disk flow, the results are discussed in Section 3, and the conclusions are presented in 
Section 4. 
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2. FORMULATION OF THE PROBLEM 

Consider an infinite disk rotating about its axis with angular velocity £2 0 (Fig. (1)). 
We take cylindrical coordinates r , 9, z with z =0 being the plane of the disk and assume 
that the fluid occupies the half space z > 0. Let p, u, v, w denote the steady-state values of 
pressure and velocity in the r , 0, z directions, respectively, in the rotating coordinate 
frame. Von Karman's exact solution of the Navier-Stokes equations for steady laminar 
rotating-disk flow takes the form 

u = r*Q 0 F(z) , v = r*Q 0 G{z) , w = (v£2 0 ) V2 H(z) , p = pv£l 0 n{z) (2.1) 

where z = z*(& 0 / v) 1/2 . The Navier-Stokes equations reduce to the following equations for 
F, G,H and P: 

F 2 -(G + 1) 2 + F'H - F" = 0 (2-2) 

2F(G + 1) + G'H - G" = 0 ^ 2 - 3 ) 

n'+HH’-H" = 0 (24) 

2F + H' = Q (25) 

where the prime denotes differentiation with respect to z. The boundary conditions are: 

F = 0,G = 0,H=0 Cz = 0) ( 2 - 6 ) 

F= 0 G = -l (2->»») (2 - 7) 

In the following analysis the variables are nondimensionalized by: 


velocity - i 

r e&0 

length 

r V "1 
^0 > 

pressure - prf Q 


V "i 

time 



1/2 


1/2 


r e n 0 


The Reynolds number is defined by 


where r* is a typical radial length scale. 

Let us write the instantaneous nondimensional velocities u , v , w and pressure p as 
u{r,d,z,t) = ^-F{z) + U{r,d,z,t) , (2.9) 

ti 


v(r,d,z,t) = -^G(z) + V(r,d,z,t ) , 
rt 


( 2 . 10 ) 


w{r,e,z,t) = — H{z)+W{r,6,z,t) , 

ri 


p(r > 8, z,t) = n(z) + P(r, 0, z y t) , 

K 


( 2 . 11 ) 

( 2 . 12 ) 


where (7, V, W are nondimensional perturbations from the steady mean flow. We substi- 
tute these expressions into the Navier-Stokes equations in the rotating frame and obtain the 
following nonlinear perturbation equations for U, V, W, and P. 


4 w+w .yl tl ,i2 + vw +w w 

I 9t R dr RdORdzRR R r dr r d6 dz 


dP_ i [ dHj_ i dHj i du d 2 u u 2 dv) 

dr R dr 2 r 2 dO 2 + r dr dz 2 r 2 r 2 dd J ’ 


(2.13) 


^L F ^OWHWF Vt r WG , + 2 (0+1)u+ UV + u W + VW w W 

at R dr R dd R dz R R R r dr r dO dz 


ldP 1 \ d 2 V 1 d 2 V 1 dV d 2 V V 2 dU 

r dd R dr 2 r 2 dd 2 r dr + dz 2 r 2 + r 2 dd 


(2.14) 


^ + L F ^ + GM HdW l dW VdW dW 

dt + R * dr + R dd + R dz + R WH +U dr + r^d +W ~& 


dP 1 \ d 2 W 1 d 2 W 1 dW d 2 W 

’ + - o -I — _ o H 1 h 


dz R dr 2 r 2 dd 2 r dr dz 2 


(2.15) 


— + — + 1— + — = 0 
r dr r dd dz 


(2.16) 


Here r is the nondimensional radius and t is the nondimensional time. We must solve 
(2.13) - (2.16) subject to the no-slip condition at the wall, whereas sufficiently far away 
from the wall we insist that the disturbances vanish. We are going to investigate the high 
frequency lower-branch instability mode for which the wavenumber varies like R' 1 at high 
Reynolds numbers as discussed in BM. 
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We introduce a disturbance into the steady motion first by perturbing the speed of 
rotation of the disk to 


£2 = i2 0 [\+Ae~ ia,r ^ , 


(2.17) 


and second by introducing a hump on the surface of the disk (Fig. 1). Here of is the dimen- 
sional frequency of the oscillations which may, for example, simulate an acoustic distur- 


bance. The shape of the hump is defined by 



( V' 2 

( * \ 

II 

* 

N 

~) Sf 

r —,9 

1 

l«0 J 

U J 


(2.18) 


where /is a function of (r* /r*) and 0; 8 and A are small scalar parameters with 5 , A « 1. 
In nondimensional form (2.17) and (2.18) become 

£2 = n 0 (l+A e~ io>,Rt ) , (2.19) 

and 

z=Sf {i' e ) • <220) 

where (0 is the nondimensional frequency defined by (0 = (0* I Since we are seeking 
solutions which are linear in 5, we can look for Fourier series solutions in the 9 direction. 
We expand f(r/ R, 9) in Fourier series as 



This will permit us to perform the analysis in the Fourier space p. We will assume that 
Reynolds number R is large and seek solutions which are linear in S and A. The perturba- 
tions U, V, W and P are written as 

U = A-^F(z)e~ i0l/Rt +5 ^{^’ z ) e ‘ Pe 

+ ASu(^,z)e- ia,Rt e ipe 

+ higher order terms ? (2.23) 

+A . 5 dL tZ y™< R ‘ e w 

+ higher order terms > (2.24) 
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1 

I 


W = A-^H( 2 )e- i0>/Rt 

R 


P = A-^-P(z)e~ ia,/Rt 

R 2 


We identify ( F,G,H,P ) as the flowfield induced by the oscillations of the disk, ( F,G,H,P ) 
as the flowfield induced by the hump, and ( u,v,w,p ) as the flowfield produced by the 
interaction between the first two flow fields. Let us assume that the factor r = r / R is 0(1). 
Substituting the above expressions (2.23) - (2.26) into the equations (2.13) - (2.16) and 
collecting terms of the same order we obtain a sequence of equations for the induced fields. 


At order A, we obtain 

(-id) + 2 F)F + HF' - 2 (G + 1 )G + F’H - F” = 0 (2.27) 

(-ico + 2 F)G + HG’ + 2(G + 1)F + G'H - G” = 0 (2.28) 

(-ico + H')H + HH' + P-H" = 0 (2.29) 

2F + H' = 0. (2.30) 

The boundary conditions at the solid boundary are 

F(0) = i7(0) = 0 (2.31) 

G( 0) = 1 . (2.32) 




6 


and the real part of /z is less than zero. For order 5 we obtain the following equations for 
F,G,H and P: 


nr* nr* n2 p 

(2 F + ipG)F + rF^r + H^--2(G+l)G + F'H = ^-Ar , 
dr dz dz 


(2F + ipG)G+rF^ + H^- + 2(G+l)F + G'H = ^- , 
dr dz ^ 


dz 1 


dH ar ixT tt dH tj,tj dP' d 2 H 
rF—— + ifiGH + H — — + H H = — — + —5- , 
dr dz dz dz 2 


„= _ dF . dH n 
2F + r — - + iBG + —— = 0 . 

dr dz 


The appropriate boundary conditions are 


(2.36) 

(2.37) 

(2.38) 

(2.39) 


F(0) = -F’(0)f(r,P) 
G(O) = -GW(r,0) 
H{ 0) = 0 


(2.40) 

(2.41) 

(2.42) 


The solution of Eq. (2.36-2.42) is: 
Fiz) = -F'iz)fir,p ) , 

Giz) = -G'iz)fir,P) , 


Hiz) = \ 2 fir,P) + r^\Fiz) + ipG(z)f(r,p) 


(2.43) 

(2.44) 

(2.45) 


and P(z ) can be obtained from Eq. (2.38). At order A5 f we obtain the equations for the 
induced disturbances £,5,iZ>andp: 


d u 


. a (9z2 TT du 


. r-'nt — L/ - i 


- (-t© + ipG + F)u-fF H A~-~ rwF’ + 2(G + l)i) = •{ -2rGG + rF| F + r ^ 


dz 


dr 


_ — - dF —dF 

+ rFF + WGFr + rH^- + rH^~ 
dz dz 


(2.46) 


~ H® + WG + F)6 - rF ^ - H - rG'w - 2 (G + l)u = f{FG + FG) + rF\ 

fa* V rfr h7 ' / 


77 -<* 3 ^ 

0+r ^, 
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(2.47) 


+fFG + ipGGr + rH— + fH— > 
dz dz 


dp d 2 w 


^- = ~^- + (-i(0 + ipG + H')w + rF^r + H^- + rF^- + ipGG+H— + H^- . (2.48) 
dz dz 2 dr dz dr dz dz 


- _ <?U _ dw n 

u + r—+ipv + r — = 0 . 

dr dz 


(2.49) 


The appropriate boundary conditions are 


u(0) = -F’(0)rf(r,P) , 


(2.50) 




(2.51) 


w( 0) = 0 . 


(2.52) 


It is easy to verify that if we write 


u = -F'(z)rf(r,p) + u x , 


(2.53) 


v = -G’(z)rf(r,p) + v 1 , 


(2.54) 


w = \2f + r^\p(z) + ipGf + w 1 , 


(2.55) 


then Eq. (2.46), (2.47) and (2.49) take the form 


• = (- 16 ) + ipG + F)u x + FF % + H + rF'w x - 2 (G + l) Dl , 

dr dz 


(2.56) 


= (- 16 ) + ipG + F)v x + rF + H + rG’w x + 2 (G + l)u x , 

dr dz 


(2.57) 


- dw , _<9u, ._ 

r^ r L + u 1 + r— i- + i)3ui = 0 , 

dz dr 


( 2 . 58 ) 


with the boundary conditions 

u 1 (r, 0) = 0 , 


(2.59) 


v 1 (r,0) = 0 , 


(2.60) 


w l (r,0) = -ipf(r,p) . 


(2.61) 


The boundary conditions at the far field are 
u[-pUi = 0 


(2.62) 


8 


v'x-pVi =0 


(2.63) 


where 

_ #(«>)- Vff(~) 2 -i'(G>+)3) 
>J ’~ 2 


(2.64) 


and the real part of p is less than zero. Equation (2.48) can be integrated to solve for 
pressure p. We observe that the above system of equations is of the parabolic type. The 
inhomogenity appears as the normal boundary condition u; 1 (r,0) = “/j3/(r,/3). Therefore, if 
we specify co f fi and /(r,/l) we can integrate these equations marching in 7 . The equations 
were solved using the two-point fourth-order compact scheme in the z direction and second- 
order upwind scheme in r direction (Spall & Malik (1989)). 

We note that the system (2.56-2.58) permits eigensolutions of the form 


V 


u Q (z) 

1>1 

■ = r 7 - 

v 0 (z) > 

w l. 

1 

Wq(z) 



r . 


(2.65) 


where yis a complex number and (uq,v 0 ,w 0 ) satisfy, 

Y- = {-i(0 + ipG + F)u 0 + jFu 0 +H^- + w 0 F'-2(G+l)v 0 , (2.66) 

dz dz 

?2 i 

— r-rp- = (-ico + i@G + F)v 0 + ~)Fv 0 + H + w 0 G' + 2(G + Dmq , (2.67) 

dz dz 


dw 0 

dz 


+ «o + yUo + iPv 0 =0 , 


( 2 . 68 ) 


with the boundary conditions u 0 (0) = v o (0) = w o (0) = 0 and u, v and w decay at infinity. 
Therefore, y is the eigenvalue of the system and determines whether the disturbance grows 
or decays algebraically in the radial direction. If in the parallel flow linearized stability 
equations (2. 16-2. 19) of BM, we replace a = a 0 / R f fi = /J 0 / R, c5 = co / R and neglect all the 
1/R 2 terms, we obtain the same equations as above (2.66-2.68) with ia replaced by y. 
Therefore, the solutions obtained for the high frequency (1/R) mode in the earlier paper are 
to be interpreted as algebraic and not as exponentially growing/decaying solutions. In 
fact, the linearized Navier-Stokes equations in a rotating coordinate system (equations 
(2.11-2.14) in BM) have solutions of the form 
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u = r 7 \u 0 (z) + 


Uiiz) 


r 2n 


v = r 7 \v 0 (z) + ^^+ 


, + ikM + . 

r 2n 


w = r r 


[ w 0 (z) u>i(z) 

I ' 9 "*’* 

i r 


p = i|i2M + M) + ... + 4S 


*t 


2 n + 1 


+ ■■ 


This solution exists for all the frequencies 0) and all azimuthal wavenumbers /3. However, 
we will see in the next section that the real part of y is positive only for frequencies (0 > 4.58 
and for a certain range of /? values, and thus the algebraic instability is not present for 
lower frequencies. 


3. RESULTS 

First we will present the values of /obtained from the eigenvalue system (2.66-2.68). 
The results are depicted in Table 1. 


Table 1 


Q) 

WMKKMIMM WM 

7 

+4.0 

+0.0 

(-.50388,15.318) 


-1.0 

(-.48892,13.176) 


-2.0 

(-.49032,10.829) 


-3.0 

(-.54533,8.089) 

+4.6 

-2.0 

(-.03848,13.484) 


-1.0 

(-.00918,15.783) 


0.0 

(.00422,17.916) 


1.0 

(.00921,19.936) 


2.0 

(.00753,21.874) 


3.0 

(.00186,23.748) 


4.0 

(-.00675,25.570) 

+7.9 

-8.0 

(-1.3621,9.593) 


-7.0 

(1.1622,14.652) 


-6.0 

(1.6793,18.261) 


-5.0 

(1.9842,21.249) 


-4.0 

(2.2082,23.919) 


-2.0 

(2.5397,28.716) 


0.0 

(2.7884,33.074) 


2.0 

(2.9893,37.156) 


4.0 

(3.1582,41.045) 
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We observe that the real part of /is positive for frequencies larger than 4.58, however, this 
occurs only for a certain range of j3 values. For (0 = 4.6, real (y) is positive in the range of /3 
= 0.0 to 3.0 and for to = 7.9, this range is from -7.0 to 40. These results are in agreement with 
the BM results reproduced in Figure 2 which shows that 1 /R mode instability is present only 
for frequencies higher than 4.58. The unstable modes computed by BM for frequencies 
lower than 4.58 do not belong to this class of instability. 

Now we present the results for the receptivity coefficient associated with the 1 JR mode 
discussed above. First we consider the case where f(r,f}) takes the form 

f(r,P) = e-’* , 

with 

* = V<y(r-r 0 ) . 

Here ^ is the location of the hump and will be taken as 7 0 =1, which is equivalent to defin- 
ing the r* based on the radial location of the hump. The parameter cr determines the spread 
of the hump in the radial direction. 

Figures 3, 4 and 5 show the downstream development of the amplitude of u lt and w l 
(Eqs. (2.56-2.58)) for the case ft) = 7.9, (i = -7 and a = 40.0. In these figures, the vertical axis 
shows the normal coordinate z and the different curves are plotted at different 7 = .6, .7, .8, 
.9, 1.0, 1.1, 1.25, 1.3, 1.4, 1.6 and 1.9. We notice that the initial form of the radial distur- 
bance velocity component u x shows two maxima, and the lower maximum disappears 
while the amplitude of the upper maximum increases continuously with 7. The 
distribution of the amplitude of the normal velocity component shows that at the edge of 
the boundary layer it approaches a constant value which increases with increasing 7. 

In Figs. 6, 7, 8, and 9 we plot the maximum amplitude of the velocity components u lt v x , 
and u>! against the coordinate x for various values of a with /? = -7 and co = 7.9. For u 1 we 
take the maximum amplitude at the outer maximum point and for u>j we take the maxi- 
mum amplitude as the amplitude at the edge of the boundary layer. We observe two differ- 
ent types of curves for a fixed a. For an x less than about .5, the amplitude pattern appears to 
be Gaussian shaped, and beyond that, the amplitude increases with increasing x. We 
identify that the first part is due to the direct influence of the wall and the second growing 
part is due to the instability. Figure 9 is the same as Fig. 7 and the curve A = 0^,(7, co)r 7 is 
plotted for a- 10, 20, 100, and 200. Here y is the eigenvalue obtained from Eq. (2.66-2.68) for 
c o = 7.9 and /3 = -7 and is equal to (1.1622, 14.652). C^P.a.co) is a measure of the receptivity 
coefficient of these modes to the excitation we considered in this paper. From Fig. 9 it is 
clear that there exists downstream an algebraic instability of the form rr and there are no 
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other eigenfunctions which exist for the system (2.56-2.58). The receptivity coefficient 
Ci(P,<7,0)) is given in Table 2 for ft) = 7.9 and P = -7 for different <7 values. 


Table 2 


<T 

Ci(P = -7, o,(o= 7.9) 

IC X I 

5 

( .00135, .00324) 

.00352 

10 

(-.00391, .00983) 

.01058 

20 

(-.04381, .01251) 

.04556 

40 

(-.11258,.01124) 

.11314 

100 

(-. 16154,.02733) 

.16384 

200 

(-.14827,. 03586) 

.15254 

1000 

(-.08145, .02089) 

.08409 

2000 

(-.05812,-01693) 

.06054 


The receptivity coefficient increases with a up to a = 100 and then decreases. Therefore, we 
conclude that there exists an optimum size for which instabilities are excited the most. In 
Fig. 10, we plot the maximum amplitude of Vi for © = 7.9, a= 100 for different P values. The 
figure shows that the waves are excited for the range of p values p = -7 to P = 20. The 
receptivity coefficients are larger for negative P values than for positive P values. 

Next, we will present the cumulative wave pattern produced downstream of the rough- 
ness site. We will consider the case ft) = 7.9 and a = 100 and the roughness element defined 
by 

f{r,e) = e-^S{d) , 
where 

x = 4o(r- 10) , 

and S(9) defines the shape in the azimuthal direction. We consider the triangular shape 
defined by 

S(6) =— (9 + a) ; -a<6< 0 

a 

= r«>-0) ; o <e<b 

O 

where a, b are positive numbers which are in the range 0 < a, b < it. If a equals b we get the 
symmetric shape; otherwise, we get the asymmetric case. We will present the results for 
the cases a = b = .05. The Fourier transform of S(0) is 
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S(0) = £s{/ 3)e ipe 


where 


S(fi) = —(a + b) 

2 it 




1__1 I l-e 
2tt p 


P = o 


j8*0 . 


Figure 11 shows the contour plot in (r,0) plane of the real part of i^(r ,0,z,t) at the height z - 
1.5 for the case co = 7.9 and < 7 = 100, a = .05 and b = .05. The roughness element is located at 
r 0 = 1 and is denoted by the point x. Near the roughness we see the effect of the roughness, 
and downstream of the roughness we observe the instability waves emerging from the 
wake. Contours of zero amplitude are highlighted by thick lines. These waves are 
inclined at angles of about -26°. The wavelength in the radial direction is about 1.0 and the 
number of waves in the azimuthal direction fi is about 6. 

In Fig. 12 we plot the maximum amplitude of IV t l for a = 4.6 for different /? and a 
values. We see that there is only weak instability for this frequency for certain values of /?. 
Figures 13 and 14 show the results for frequencies a = 4.0 and 0.0. For these frequencies, 
all the disturbances decay downstream. All these conclusions agree with our previous 
calculations for yand the results of BM that the 1 /R instability is present only for <u > 4.58. 


4. CONCLUDING REMARKS 

The receptivity and non-parallel stability of travelling disturbances of long 
wavelength (on the order of Reynolds number R) are investigated. The non-parallel 
effects are included in a self-consistent manner. It is shown that eigensolutions exist 
which grow or decay algebraically in the radial direction. It is established that the 
solutions obtained for the high frequency (1 JR) mode in the earlier paper of BM are to be 
interpreted as algebraic and not as exponentially growing/decaying solutions. These 
algebraic-type solutions exist for all frequencies (O and all spanwise wavenumbers /J. 
However, growing disturbances exist only for frequencies larger than 4.58 times the disk 
rotational frequency and for a certain range of /3 values which is in agreement with the 
results of BM. 

The receptivity calculations show that there exists an optimum size of roughness for 
which instabilities are excited the most. For example, the receptivity coefficient Ci for 
frequency (O - 7.9 with spanwise wavenumber j3 = -7 is .01, .16 and .084 when a = 10, 100 and 
1000, respectively. Gaussian shaped roughness elements with a = 100 correspond to a 
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roughness element with a spread of about . 1 r 0 in the radial direction, where r 0 * is the 
location of the roughness. This shows that very narrow or broad roughness elements will 
not excite these modes efficiently. It is also observed that the receptivity coefficients are 
larger for negative /J values than for positive fi values. Hence, if these modes are excited 
they will align along negative angles. 

The cumulative wave pattern produced from an isolated roughness element shows that 
the waves are inclined at about -26° and the wavelength in the radial direction is about 1.0. 
In dimensional variables, the wavelength 1.0 equals to /q where the is the location of the 
element. In some experiments performed by Faller and Kaylor (1966), in Ekman layer 
and rotating disk boundary layer, the waves oriented at negative angles have been 
observed at finite Reynolds numbers. These waves were rapidly moving and had long 
wavelengths compared to cross-flow vortices. Even though there seems to be some 
similarities between these computed algebraic modes and the experimental observations, 
more measurements have to be done to draw any definite conclusions. 

Another remaining question is how the lower branch exponential-type solutions join 
with this algebraic mode. For example, consider the neutral stability curve for frequency 
(O = 6.0 shown in Fig. 2. The lower branch neutral curve, which is asymptote to the neutral 
curve for stationary disturbances, is an exponential-type solution. We see that this curve 
turns around at a Reynolds number of 3x10^ and becomes a 1 JR mode. However, we 
discussed that this 1 JR mode is actually an algebraic-type mode. Therefore, why does the 
lower branch neutral curve, which is asymptote to the neutral curve for stationary 
disturbances, turn around, and what is its subsequent behavior, are the remaining 
questions. 
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Fig. 1. Schematic of rotating-disk flow. 
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Fig. 2. Lower branch of neutral stability curve for rotating-disk flow: (a^-plane using 
log-log scale. 
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Fig. 10. 



26 


Pig. 11. Contour plots of V* velocity in (r\0) plane at 
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Fig. 12. Maximum amplitude of V\ velocity for the case (0 = 4.6 at different /3 and o 
values. 
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Fig. 14. Maximum amplitude of 
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